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Abstract 

Let (M, g) be an analytic, compact, Riemannian manifold with bound¬ 
ary, of dimension n > 2. We study a class of generalized Radon trans¬ 
forms, integrating over a family of hypersurfaces embedded in M, satis¬ 
fying the Bolker condition [^. Using analytic microlocal analysis, we 
prove a microlocal regularity theorem for generalized Radon transforms 
on analytic manifolds defined on an analytic family of hypersurfaces. We 
then show injectivity and stability for an open, dense subset of smooth 
generalized Radon transforms satisfying the Bolker condition, including 
the analytic ones. 


1 Introduction 

Let (M, g) be an analytic, compact, Riemannian manifold with boundary, of 
dimension n > 2, with volume form denoted by dVol. Let E be a family of em¬ 
bedded hypersurfaces. A generalized Radon transform takes each / £ C°°{M) 
to the set of its integrals over the hypersurfaces of E, with respect to the sur¬ 
face measure induced by the volume form. Often E is itself a smooth manifold; 
for example, the Euclidean Radon transform is defined over the set of affine 
hypersurfaces in R", which has M x 5'"“^ as a double cover. Such transforms 
are found in applications to many other fields, including harmonic analysis, 
scattering theory, seismology and medical imaging. 

The main questions regarding these transforms include determining condi¬ 
tions under which they are injective, finding when the transform has a stable 
inversion, and characterizing the range. We concentrate on the first two ques¬ 
tions here. These questions are also important in the partial data case, where 
integrals are known only for a subset of E. In the case of the Euclidean Radon 
transform, we refer the reader to [M m [HI HU and references therein for the 
resolution of these problems, and their generalization to the Radon transform 
over symmetric spaces and other related contexts. 

The geometric data of a generalized Radon transform can be encoded by an 
incidence relation between points on M and the hypersurfaces in E that contain 
them. Let A C M x E be this relation, i.e., the set of ordered pairs (x, cr) such 
that X € a. One says A is a double fibration when it is a smooth, embedded 


1 


submanifold of M x S such that both canonical projections are smooth and 
their restrictions to A form a fiber bundle over M and S, respectively [lOj . 

Guillemin and Sternberg [HIIlIIl] showed that given a Radon transform 
R defined by a double fibration A, both R and its adjoint R* (often called 
the generalized backprojection by analogy with the Euclidean case) are Fourier 
integral operators, and the canonical relation of R is the conormal bundle N*A 
of the incidence relation. If in addition A satisfies the Bolker condition, which 
says that the induced projection tt* : N*A —>■ r*E is an embedding, then the 
normal operator R*R is an elliptic pseudodifferential operator, which yields 
invertibility up to smoothing error. 

A stronger result than invertibility is a Helgason-type support theorem, by 
analogy with that of the Euclidean Radon transform m- Such a support the¬ 
orem implies that if / is a priori of compact support, and Rf = 0 for all 
hypersurfaces intersecting the support, then / = 0. In the analytic category, 
there has been much work in this area (for example, [31 01 [331 [33]) using an¬ 
alytic microlocal analysis and the Bolker condition to prove support theorems 
for analytic generalized Radon transforms (i.e., with M, S analytic manifolds) 
with nonvanishing analytic weight. On the other hand, in n = 2 there is a coun¬ 
terexample in the smooth category due to Boman [3] of a function supported in 
the disk such that some weighted Radon transform over lines vanishes. For the 
weighted X-ray transform over curves, there is an analogous Bolker condition 
m and support theorems are known for a class of such transforms in n > 3 
[311133] including the geodesic ray transform on an analytic, simple manifold 
over functions [20] and over symmetric tensor fields [3Tj . 

Our first result considers such analytic generalized Radon transforms and 
shows their analytic microlocal regularity. This builds upon similar results for 
the weighted X-ray transform over a generic class of curves BIS] and in particu¬ 
lar, geodesics [371 [28]. To avoid complications at the boundary of M, we embed 
it isometrically in a slightly larger, open analytic manifold Mi. In section 2, we 
show how to extend the definition of R^ to a transform on Mi in a stable way. 
We show: 

Theorem 1. Let R^, he an analytic generalized Radon transform satisfying the 
Bolker eondition, with w an analytic, nonvanishing weight. Let f G £'{Mi) be 
sueh that Rwfio') = 0 in a neighborhood of some hypersurface (Tq G E. Then the 
analytic wavefront set WFa(/) does not intersect the conormal bundle N*ao. 

The main tool is a complex stationary phase lemma of Sjdstrand [261 Theo¬ 
rem 2.8, 2.10 ff.] and related techniques, which suffice in lieu of a hypothetical 
analytic calculus of Fourier integral operators. The proof of this theorem is 
given is section 3. 

Theorem [T] also implies a local support theorem [231 Prop 2.3] and in par¬ 
ticular the injectivity of : L^{M) —>■ In fact, the proof of 

the theorem shows that it suffices for Ru,f{(j) to be analytic in a neighborhood 
of (To- One then obtains a unique continuation result of the following type: if 
/ S £'{Mi) is analytic on one side of cto , and R.u,f{a) is analytic in a neighbor¬ 
hood of (To, then there is a neighborhood of (Tq on which / is analytic. 
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Our second result is a stability estimate for a generic class of smooth gener¬ 
alized Radon transforms satisfying the Bolker condition. We restrict ourselves 
to those generalized Radon transforms studied by Beylkin [T], which have E 
parametrized globally by the level sets of a smooth defining function ip satisfy¬ 
ing some conditions to be made explicit later. 

Theorem 2. Let {M,g) be an analytic Riemannian manifold with boundary. 
Take Rw : L^{M) —>■ to be an injective generalized Radon trans¬ 

form defined by ip with weight w, satisfying the Bolker condition. Then there 
exists K ^ n and a neighborhood of {ip, w) € such that the generalized 
Radon transform Rw defined on (M, g) by a defining function and weight in this 
neighborhood is injective and for all f G L'^(M) there exists C > 0 such that 

This follows from an analysis of the symbol of the normal operator . As 

mentioned above, under the Bolker condition it is an elliptic pseudodifferential 
operator. We show that perturbing the defining function and weight slightly in 
perturbs the operator slightly, preserving the stability estimate. 

While we work entirely on an analytic Riemannian manifold {M, g) and 
do not perturb the metric in this result, we use the metric only to provide a 
convenient choice of surface measure, and to ensure the existence of a dense set 
of injective, generalized Radon transforms. We may then conclude: 

Corollary 1. On each analytic, compact Riemannian manifold with boundary, 
there is a generic set of generalized Radon transforms satisfying the Bolker con¬ 
dition that are both injective and stable. 

We defer the proof of Theorem [5] to section 4. 

Ackno-wledgements. The first author is partly supported by NSF Grant DMS- 
1301646. 

2 Generalized Radon Transforms 

In this section we fix notation and establish some basic facts about the gen¬ 
eralized Radon transform, including a statement of the Bolker condition. For 
concreteness we consider the space of hypersurfaces S as parameterized by a 
defining function, following Beylkin [1], though we only consider oriented hy¬ 
persurfaces. To avoid difficulties occuring at the boundary of M, we assume M 
is isometrically embedded in a slightly larger open manifold Mi, whose metric 
we also refer to by g. If we are considering the analytic category of Radon 
transforms, we will also assume Mi is analytic. In the sequel, we will always 
consider L^(M) to be functions on Mi, extended by zero. 

Definition. Let ip € C°°{Mi x (R” \ 0)). p is a defining function when it 
satisfies the following conditions: 
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2 . 

3. 


ip{x, 9) is positive homogeneous of degree one in the fiber variable, 
if is non-degenerate in the sense that dx<pix,9) ^ 0. 

The mixed Hessian of tp is strictly positive, i.e., 

d'^ip 


det 


dx'^dOi 


> 0 . 


The level sets of ip will be denoted by 

Hs,e = {x G Ml : (p{x,6) = s}. 

Note that by homogeneity, Hg^e = H\s,\e for A > 0. Therefore we can 
consider E as globally parameterized by {s,6) € M x Often we will also 

implicitly consider as a function on Mi x 

The third condition imposed on a defining function is a local form of Bolker’s 
condition. This allows us to locally identify {x, 9) G Mi x 5'”“^ with the covector 
dxip{x,9)/\dxip{x,9)\g G S*Mi. We will assume in addition a stronger, global 
Bolker condition. 


Definition. A defining function ip satisfies the global Bolker condition if for 
each 9 G the map x i—>■ dgip{x,9) is injective, and for each x G M, the 

map 9 I—>■ dxip{x, 9) is surjective. 

The first condition is roughly analogous to the “no conjugate points” condi¬ 
tion assumed by 0110] for similar results regarding the geodesic ray transform, 
and the second ensures that every singularity is observable from some hypersur¬ 
face in E. Note that generalized Radon transforms defined by a double fibration 
satisfying the Bolker condition as stated by Guillemin et. al. also satisfy this 
Bolker condition, see [23l Lemma 3.5]. 

Consider / G C°°{M). We extend it by zero to a function on Mi which 
we also denote by /. Let the generalized Radon transform i?u, determined by 
(M, g, ip, w) be defined by 

Rwf{s,9)= ( w{x,9)f{x)dps,e, 

JH,,g 

where w G C°°{Mi x 5'"“^) is a smooth, nonvanishing weight and dps.e is the 
volume form on H^ g induced by dVol. There exists a smooth, nonvanishing 
function J{x,9) such that 


dpLsfi{x) Ads = J{x, 9) dVol. 

We calculate the adjoint of Rw in L^{M,dYo\) to be 

[ [ {Rwf)gdsd9 = ( f f w{x,9)f{x)g{s,9)dp.s,edsd9 

dS"--i Jr Js^-^JRJHsg 


I Ml 


g{ip{x, 9), 9)w{x, 9)J{x, 9)f{x) dVol d9 
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Therefore 



w{x, 9)J{x^ 9)g{(p(x, 9), 9) d9. 


This is simply a generalized backprojection with weight wJ. 

3 Microlocal regularity 

In this section, we take {Mi , g) to be an analytic Riemannian manifold, ip to 
be an analytic defining function, and w to be an analytic nowhere vanishing 
weight. Given / £ £'{Mi), we are interested in the microlocal analyticity of 
/ given that of Rwf- (We extend Rw to £'{Mi) by duality.) We will use the 
following definition of the analytic wavefront set, following Sjdstrand. There are 
alternative approaches to analytic wavefront set by Sato, Kawai, Kashiwara [25) 
and also Bros and lagolnitzer [B], which were shown to be equivalent by Bony 


0 . 


Definition ([SQ Def. 6.1]). Let {xo,^o) € T*]R” \ 0 and let ij}{x,y,^) be an 
analytic function defined in a neighborhood U of (xq, xq, ^o) £ such that 

1. For all {x,x,ff) £ U (i.e., x = y), we have 


tp{x, a;, ^) = 0 and 


2. There exists G > 0 sueh that for all {x,y,^) £ U, we have 

lm'il;{x,y,() > C\x - yp. 


Let a{x,y,^) be an elliptie elassical analytic symbol defined on U, see, e.g., 1261 
Theorem 1.5]. 

We say u £ 2?'(IR”) is analytic microlocally near (a;o,Co) if there exists a 
cut-off funetion x G G)]”(M") with xi^o) = 1 such that 



for some C > 0, uniformly in a conic neighborhood of (xo,^o)- 

The analytic wavefront set is the closed eonie set WFa{u) C T*K”\0, which 
is the complement of the set of covectors near which u is microlocally analytic. 

We note that this definition is microlocal and invariantly defined, and there¬ 
fore can be extended to distributions on analytic manifolds (see [II Theorem 
8.5.1] and the remarks following). In this case, for u £ VfMi), WF^(u) is a 
closed conic subset of T*Mi \ 0. 

Recall that since the mixed Hessian of ip is strictly positive, we may locally 
identify {x,9) £ Mi x with the unit covector dxp{x,9)/\dxip{x,9)\g £ 

S*Mi. Fix a covector (xq, 9q) £ T*Mi \ 0 with sq = ^o)- From now on we 

will work in a small conic neighborhood of this covector. 
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Proposition 1. If R^f{s,9) = 0 for (s,0) in a neighborhood of (sqj^o); then 
{xo,d:„ip{xo,eo)) ^WF^(/). 

Proof. Let us fix a coordinate system. We already have local coordinates (x, 9) 
on T*M \ 0. Without loss of generality we can take sq = 0 and |0o| = 1- To 
simplify the coordinates on S, we perform a stereographic projection onto the 
tangent plane of the sphere at 9o, which is an analytic diffeomorphism mapping 
a neighborhood of 9q G 5'"“^ to a neighborhood of the origin in We 

refer to the coordinates on this tangent plane by and pass to a perhaps 
smaller neighborhood of E with |s| < 2e and |^| < S, with e, d > 0 being small 
parameters. 

Much of the complexity of analytic microlocal calculus is due to the dif¬ 
ficulty of localizing in the analytic category, as there are no suitable cut-off 
functions. Instead one often uses a sequence of quasianalytic cut-off functions 
Xn G C“(K), depending on e, for whose construction we refer to [HISD]. We 
will only use the following properties of this sequence: 

1. suppxAT C {-2e,2e) and xjv(-e,e) = 1. 

2. For all G N and k < N, the estimate 


<iCN) 


holds for a constant C > 0 independent of N. 

By assumption Rujf{s,f) = 0 for |s| < 2e and |^| < S. Let A ^ 1 be a large 
parameter, to be fixed later. This implies that 


0 = / e*^®XAr(s) [ w{x,^)f{x)dns,ids. (1) 

J 

Recall that f are analytic coordinates for the neighborhood of 9q in 5'”“^ that 
we are concerned with, and so here and in the sequel we write for brevity, e.g., 
w{x,f) = w{x,9{f)) and d^s,^ = dMs,e(5)- 

It follows from Beylkin’s construction that 

Ads = J(x, dVol 

where J{x,f,) is an analytic, nonvanishing Jacobian and dVol is the volume form 
on Ml associated to the metric. Hence (HJ reduces to the oscillating integral 


of^v{^,i)aj^{x,f)f{x)dYo\ = 0. 


( 2 ) 


Here a]v(x,f) is a sequence of classical analytic symbols on the same neighbor¬ 
hood of (xo,0) G Ml X The coordinates on x and f are real-analytic, 

and so we may extend their domain of definition slightly by analytic continu¬ 
ation to a Grauert tube of a small neighborhood of iLo,o C Mi (for x) and a 
small neighborhood of the origin in for f. This continuation in principle 
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depends on the choice of analytic coordinates, but as the analytic wavefront 
set is invariantly defined the final result does not depend on this choice. We 
choose a perhaps smaller 6 such that G : |^| < 5/2'\ is contained in this 
neighborhood. We denote the local complex coordinate patch of a:o as C/ C C". 

Let y G U and 77 G C"“^, with \r]\ < <5/2. Let p{^) = 1 when |^| < S and 
zero otherwise. Then we multiply (I2|) by 

p(^-77)exp , 

and integrate with respect to The resulting integral is of the form 

jj r])f{x) dVo\{x) = 0. (3) 

Here bN is a sequence of classical analytic symbols defined on a complex neigh¬ 
borhood of Hq q X {0} X {0} and 3> is the augmented phase function given by 

To estimate the left-hand side of ([3]), we intend to use the method of com¬ 
plex stationary phase. Therefore, we are interested in the critical points of the 
function f $(x, y, 77 ). Note that 

y, v) = -v) + 0 - d^^iy, 0- 

There are clearly real critical points ^ when ^ = 77 and x = y. These crit¬ 
ical points are non-degenerate, and therefore induce complex critical points 
ic{x, y,r]) =r] + i{y - x) + 0{6). 

Consider the situation when y = 0. Then for x ^ 0, the only real critical 
points are where d^(p{x,^) = d^ip{y,^). However, this cannot happen by the 
global Bolker condition that we imposed on the defining function. By non¬ 
degeneracy again we see there are no real or complex critical points other than 
^cix, y, rf) for (x, y, 77 ) where | 7 /| < 5 and |C - ??| < S. 

Now we apply the complex stationary phase lemma [26l Theorem 2.8, 2.10] 
to As a preparatory step divide the integral into two regions; one over the 
region 

1 + = {(a:, y, v) --Ix-yl^ (5/Co, IC - ?7l < <5} 

and one over the region 

I- = {ix,y,tv) ■\x-y\> S/Co,\^-v\ < <5}- 

Here Cq > 0 is a constant chosen so that the critical points ^c(a;, y, 77 ) lie within 
and none lie in 

In I- , we may define the usual operator L such that via 



This is well-defined as there are no critical points in so we may repeatedly 
integrate by parts: 


/ d^'^bNfdYold^ 

— 

f {L^d^‘^)bNfdVo\d^ 




J 



< 

f e^^'^iL*f[bNf]dVold^ 



J I- 



N 




The terms Bk are boundary terms that decay exponentially, due to the fact that 
Imd) > 0 for 1^ —? 7 | = 0((5). As for the integral on the right-hand side, we recall 
that b]s[ is defined by 


bN{x,C,r]) = p{C - r])xN{<f{x,0)w{x,^)J{x,^). (4) 

The worst possible growth of {L*)^b]sf in terms of N occurs when all derivatives 
are applied to ^)), and in this case we may apply the estimate 


di^hNis) <iCN) 


\N 


which follows from the construction of the sequence of quasianalytic cut-off 
functions. Therefore, 


dx 


' I- 


= o 




(5) 


As for the integral over /+, the cut-off functions Xn{‘-p{x,^)) are all equal 
to one. Therefore the amplitude on /+ does not depend on iV; we remove this 
dependence and refer to the amplitude restricted to this region as b. We know 
all of the critical points of ^ i—>■ and can therefore apply the complex stationary 
phase lemma. This yields an estimate of the form 

J e^^’^bfd\o\di = C\-'^/‘^j B f d\o\ +O ({CN /X)^ + Ne-^/^ 

Here 'tp{x,y,r]) = <^{x,y,^c{x,y,r]),r]) and B{x,y,y) = b{x,^c{x,y,y),y). We 
may now fix N such that N < (A/Ce) < TV -|- 1 to ensure the error is exponen¬ 
tially small. 

giAV-^/^ol ^ o(e-VC)^ 


/' 


Now B{x,y,r]) is an elliptic analytic symbol near {x,y,r]) = 0 and tp(x,y,r]) 
is a non-degenerate phase function. To show this implies (a:o,6>o) ^ WF^(/), 
we check the details of the characterization of the analytic wave front set given 
above. Recall 


i^ix, y,y) = - |^c(a:, y, y) - y\^ -f p{x, ^c{x, y,y)) - f{y, ^dx, y, y)). 












Note that ^c{x, x,r/) = rj for x real, and therefore 'fpi.x, x, 77 ) = 0. In addition 

d^ili{x, X, vj) = d^ip[x, vj) = -dy-tpix, X, rj). 

By the global Bolker condition we can make a change of variables 77 ' so that 
77 ' = dx^p{x,ri). Finally, it is clear that Im 7 /;(x, 7 /,.^) > C\x — j/p for x,y real. 
Therefore, (xq, ^o)) ^ WFa(/). □ 

Theorem [ 1 ] follows from applying the proposition to all conormals of a fixed 
hypersurface cto- 

Remark. From the proof we see that it suffices for Ryjf{a) to be analytic in a 
neighborhood of (sq) do). After microlocalization, the right-hand side of ([2|) will 
be instead of zero, but this poses no problem. 


4 Stability 

We now return to generalized Radon transforms with smooth defining function 
175 : Ml X 5'"“^ and smooth, nonvanishing weight w : Mi x The object of 

interest in this section is the normal operator Nw = R^Rw It is known that 
the global Bolker condition implies Nw is a pseudodifferential operator [131 Prop 
8.2]. However, we require more detailed knowledge of the symbol of Nm for the 
kind of stability estimates we prove later. 

First we obtain a representation of the Schwartz kernel of Rw 

Lemma 1. The Schwartz kernel G x 5'"“^ x Mi) of Rw is 

KR„{s,0,y) = (27r)"M(s - ip{y,0))w{y,0)J{y,0) 

where J{y,0) is the smooth, nonvanishing function such that 

dfj.s,e{y) /\ds = J{y,0)dVo\{y). 

Proof. We perform a partial Fourier transform of Rwf{s,0) in the s variable, 
taking s' to be the dual variable of s. The change of variables then yields 


FsRn,f{s',0) 




w{y, 0 )f{y) dfj,s,g ds 



e— 0 ) j(y, 0 )f(^y) cNo\{y). 


Therefore 


R„/(s, 0) = (2^)-i [ [ 0)J{y, 0)f{y) dVol(y) ds' 

Ml 


KR,„{s,0,y)f{y)dVo\{y). 


□ 


Ml 
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Similarly, the kernel of the generalized backprojection R'^ is 

5{{p{x,9) — s)w{x,0)J{x,0). 

From this we see that the kernel of is 

Kn^ = (27r)-i Jj 9)1 {x, 9)w{y, 9)J{y, 9) ds' d9. (6) 


We can now use this representation to find the principal symbol of the normal 
operator N^. 


Lemma 2. The principal symbol of is 




(27r 


+W{x,x,-£,/\i\) 


where W is the auxiliary function 


W{x,y,9) = w{x,9)J{x,9)w{y,9)J{y,9). 

Proof. Beginning from m, we split the integration over R into {s' > 0 } and 
{s' < 0}. Using the positive homogeneity of the defining function, we rewrite 
the integral as 

Kn„= f [ 

Jo 

dS"-i Jo 

Here and are the Schwartz kernels of the operators Nff and N~ 
respectively, so that Nyj = iV+ + N~. We work with each term separately. 
Let f = s'9 be polar coordinates for R”. This change of variables is justified 
when the kernel is applied to a test function in using the proof of [THl 

Theorem 7.8.2] it can be shown that it is justified for the kernel itself. Then we 
obtain 

" L |j) l^r”” 

By the global Bolker condition, d^ip{x,^) = d^(p{y,f,) implies x = y. A sta¬ 
tionary phase argument implies that K'^ is a smooth function away from the 
diagonal of Mi x Mi. 

Fix xq € Ml. There exists a neighborhood U of xq on which we have normal 
coordinates, which we refer to again with (a:®), such that a:(a;o) = 0. We then 
use (a;®,y®) as coordinates on U x U, with x® = y®. We consider the localized 
kernel 

xKti^X = J ^x, y, ili ^ 
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In these local coordinates, we can expand the phase function near the diagonal 
x = y. 


0 - ip{y, = {x-y)- / da:ip{x + t{y -x),^) dt 
Jo 


Define the map 


^'{x,y,^)= d^ip{x + t{y-x),^)dt. 


Jo 


Near the diagonal, this map is smooth, and 
'di' 


det 


di 


{x, a;, 0 = det id^a:^{x, 0) = C) > 0- 


and so the map {x,y,^) i—>• {x,y,^') is a diffeomorphism of a neighborhood of 
the diagonal onto another neighborhood of the diagonal. It may be necessary to 
shrink the support of x slightly for the change of coordinates to be well-defined. 
In these variables, 

= ix-y)- 

Both sides are positive homogeneous of degree one, which implies that |^'| = 
c(a;,2/)|^| with c{x,y) a strictly positive, smooth function defined near the diag¬ 
onal. Clearly c{x,x) = |cIa:V5(a;,^)|. This reduces the cut-off kernel of to an 
honest pseudodifferential operator 

xKm^X = ^a;, y, x(a^)x(2/)ICT“”c(a:, y)”"^ 

To evaluate the principal symbol of , we restrict the amplitude to the 
diagonal x = y. The principal symbol of is the sum of those for iV+ and 
N-. □ 
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This reconfirms that iV^, is an elliptic pseudodifferential operator of order 
1 — n, provided the weight is nonvanishing and the global Bolker condition is 
satisfied. 

We now consider the stability of reconstructing / G L^(Mi) from i?™/, using 
the analysis of the normal operator = R^jR^f from the previous section. 
The basic estimate follows from elliptic regularity; the stability estimate follows 
from [5ni Prop V.3.1]. 

Lemma 3. Let w G C°^{M x 5”“^) be a nonvanishing weight and let ip G 
C°°{M X (H." \ 0)) be a defining function. Then for all f G L'^(M) and s > 0 
there exists C > 0 and Cg > 0 depending on s such that 

||/|| l 2 ( m ) < C\\Nu,f\\H”-^{Mi) + Csllflln-^- 

If, in addition, ■ Lf{M) —> is injective, then we have a stability 

estimate with a loss of n — 1 derivatives, 

||/||l2(M) < C\\Ny,f\\Hn-l(Mi) 
with a different constant C > 0. 
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In particular, by Proposition 1, the latter stability estimate holds when the 
geometric data (i.e., M,ip and w) are analytic. This also follows directly from 
[3]. Our main contribution is to extend this stability estimate by perturbation 
to a generic set of smooth geometric data. We begin by using the standard pseu¬ 
dodifferential calculus to show that the normal operator depends continuously 
on finitely many derivatives of the data. 

Lemma 4. Let [Mi,g) he an open Riemannian manifold with an embedded 
compact manifold M with boundary. Let ipi,ip 2 be two defining funetions and 
wi,W 2 be two nonvanishing weights. Let Ni = and N 2 = R^^Rw 2 - 

There exists a K ^ n such that if 

WTI - T2\\c^{MixS”'-^)^ ll'W^l - W^2 ||c«(MixS"-i) < ^ < 1, 

then there exists C > 0 depending a priori on the {Mi x 5”“^) norm of (pi 
and wi such that 


\\{Ni-N2)f\\H^-^M,)<C6\\f\\L^M,)- 

Proof. We have seen in the previous lemmas that Ni and N 2 are both elliptic 
pseudodifferential operators with symbols depending on ipi , ip 2 and wi , W 2 re¬ 
spectively. Let K be an arbitrary, large natural number to be fixed later. If 
the defining functions and weights are d-close in {Mi), then it follows from 
Lemma 2 that the amplitudes are 0{S) in (7^“^(Mi). By the continuity of 
pseudodifferential operators m Theorem 18.3.II and ff.] the operator norm of 
is bounded by a constant multiplied by some C'°°(Mi)-seminorm of 
the difference of the amplitudes. Take K large enough so that 

llL2(Mi)->ff'*-i(Mi) = 0(^)- 

The lemma follows from adding the positive and negative parts of the estimate 
together. Notice that K, the necessary number of derivatives, does not depend 
on the defining functions or the weights themselves. □ 

It is of interest to determine the minimal regularity necessary for the per¬ 
turbation result of the previous lemma. For the related geodesic ray transform, 
it is known that the geometric data need only be i5-close in [9]. One would 
then expect the above to hold for data d-close in C". 

We now able to prove with our main result. Theorem [31 

Proof of Theorem\^ Recall is an injective generalized Radon transform. 
Lemma 3 yields the following stability estimate: 

||/||l2(Mi) < Oilli?* i?^,/||//n-l(Ml)■ 

Then Lemma 4 allows us to perturb this estimate using 
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Therefore, 


For S < min{(2CiC2)“^, 1/2}, the second term on the right-hand side may be 
absorbed into the left. The resulting stability estimate for the perturbed normal 
operator implies injectivity of Rw. □ 
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